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The density profile p{z) and Gibbs adsorption T of a near-critical fluid confined between two 
i— i, identical planar walls is studied by means of Monte Carlo simulation and by density functional 

• theory for a Lennard- Jones fluid. By reducing the strength of wall-fluid interactions relative to 

^ ' fluid-fluid interactions we observe a crossover from behaviour characteristic of the normal surface 

S universality class, strong critical adsorption, to behaviour characteristic of a 'neutral' wall. The 
crossover is reminiscent of that which occurs near the ordinary surface transition in fsing films 
! 1 * ' subject to vanishing surface fields. For the 'neutral' wall p(z), away from the walls, is almost 

constant throughout the slit capillary and gives rise to a V that is constant along the critical isochore. 
ryj ' The same 'neutral' wall yields a line of capillary coexistence that is almost identical to the bulk 

! t . coexistence line. In the crossover regime we observe features in p{z) similar to those found in the 

magnetisation profile of the critical Ising film subject to weak surface fields, namely two smooth 
maxima, located away from the walls, which merge into a single maximum at midpoint as the 
strength of the wall-fluid interaction is reduced or as the distance between walls is decreased. We 
discuss similarities and differences between the surface critical behaviour of fluids and of Ising 
C | magnets. 

, Q, l PACS numbers: 64.60.Fr, 05.70.Jk, 68.35.Rh, 68.15.+e 

> ■ 

00 1 I. INTRODUCTION 

<n : 

When a fluid, a fluid mixture, a binary alloy or an Ising magnet that is close to a critical point in bulk is brought 
i to a surface its properties are altered significantly. In particular the bulk universality class is divided into several 

£f) ' surface universality classes which depend on whether the tendency to order at the surface is increased or decreased 
w.r.t that in bulk [1,2]. For the case of Ising magnets the phenomenology is well- understood: the critical behavior 
of the semi-infinite system is governed by a surface scaling field c which describes the enhancement of interactions in 
the surface layer, c > corresponds to a reduced tendency to order in the surface and c = oo defines the fixed point 
of the renormalization group transformation of the ordinary transition. Similarly c = — oo defines the fixed point 
for the extraordinary transition and c — c* corresponds to the special transition. Such a classification pertains to 
vanishing (external) surface field hi = 0, a situation which is easily realized for the Ising magnet but not for a fluid; 
the containing walls must always exert a non-zero surface field of some type. Usually the wall or substrate is assumed 
to exert an effective potential on the fluid that is infinitely repulsive at short distances and strongly attractive at 
large distances. Such potentials give rise to pronounced peaks in the density profile of the fluid near the substrate 

• i— i , leading to strongly positive adsorption and in magnetic language this corresponds to hi ^> 0. On the other hand, 
the presence of a substrate should decrease the net fluid-fluid attraction between an atom and its nearest neighbors 
below the corresponding bulk value so that c > 0. Thus fluids are expected to lie in the universality class of the 
normal transition with fixed point hi = oo and c — oo. Since the shapes of the order parameter profiles should be the 
same whether the symmetry is broken by the external surface field, hi > 0, or spontaneously, c < 0, the normal and 
extraordinary transitions should be equivalent [3,4], provided the symmetry breaking fields are sufficiently (infinitely) 
strong. 

In recent years much effort has been focused on understanding the crossover from the normal to the ordinary 
transition, i.e. the regime of weak surface fields. The length scale associated with the surface scaling field hi is 

li ~ h 1 "/ Al 5 where v is the bulk correlation length exponent and A° rd is the surface gap exponent. For a semi- 
infinite Ising-like system Ritschel and Czerner [5] showed that for small hi the order parameter profile can take a 
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non-monotonic form, increasing with distance z from the surface for z < l\ before decreasing for z>l\. This is quite 
different from the case of strong field adsorption where the profile decays monotonically, as z~P/ v ', at the critical point. 
Such behaviour is not predicted by mean-field treatments of Ising systems but was confirmed in several subsequent 
studies [6-9]; it is also found in exact results for the d = 2 Ising model [10-12]. Physically l\ corresponds to the 
distance from the surface over which the order parameter responds linearly with field hi [6]. 

The existence of the critical length scale l\ has important consequence for confined systems near bulk criticality. 
There are three relevant lengths: h, the distance D between the two confining surfaces which both exert fields hi 
and the bulk correlation length £(,. Sufficiently close to bulk criticality one can investigate the regime where D is 
comparable to l\ with D,l\ <C In a recent study [13] of the critical d — 2 Ising film order parameter profiles with 
two symmetric maxima at z ~ l\ and z ~ D — l\ were obtained for weak fields hi. For even weaker fields when l\ <~ D 
the maxima merge into a single one at midpoint, D/2. That study provides some of the motivation for the present 
one where we investigate the density profiles and the Gibbs adsorption for a simple (Lennard-Jones) fluid at bulk 
criticality that is confined by two planar walls separated by a distance D. We enquire whether phenomena analogous 
to these found for Ising systems subject to weak surface fields are found in fluids. 

As implied earlier, it is not at all obvious that one can achieve the analogue for a fluid of surface magnetic field 
hi — > 0. Thus it is not clear, a priori, that a physical length scale l\ should exist. On the other hand, one knows 
that for a purely repulsive wall complete drying with adsorption r — > — oo (wetting by gas) occurs below the critical 
temperature T c whereas for a strongly attractive wall complete wetting by liquid occurs, r — > oo, so that one might 
imagine a scenario whereby for weak wall-fluid attraction the adsorption T is close to zero. Moreover, if T w for 
all states on the critical isochore, p = p c , T > T c , one might argue that the fluid is exhibiting one of the features of 
the Ising magnet subject to a vanishing h\. In other words, by suitably tuning the wall-fluid potential can one mimic 
the situation hi — 0, i.e. the ordinary transition? Of course, finding a wall-fluid potential that corresponds precisely 
to hi — 0, i.e., one that generates precisely the same features in the profiles and, therefore, the adsorption for all 
T, seems very unlikely - given the lack of Ising symmetry in any 'real' fluid. However, if one does find a wall-fluid 
potential that is 'neutral', giving rise to T ~ 0, then one might also expect to see phenomena associated with small 
hi, i.e. the mesoscopic length scale h. The existence of such phenomena would imply that the ordinary transition is 
not unique to systems with Ising symmetry and would lend support to the general idea of fluid-magnet universality. 

We are not aware of any published proof that a fluid system cannot undergo the ordinary transition [14]. The only 
rigorous result for fluids [15] is for a planar hard wall where an exact sum rule analysis shows that the surface critical 
behaviour lies in the universality class of the normal transition and the surface critical exponents can be expressed in 
terms of bulk critical exponents, as predicted originally in Ref. [3]. 

What is the experimental situation? To the best of our knowledge there is no experiment that probes critical 
adsorption in the weak surface field regime for a one-component fluid [16]. For binary fluid mixtures there is strong 
evidence from (i) reflectivity studies of critical carbon disulphide-nitromethane mixtures adsorbed at a silanated 
borosilicate glass surface [17] and (ii) ellipsomctric studies of an homologous series of critical n alkanc - methyl 
formate mixtures [18] that the weak adsorption, hi small limit, is achieved. In the second of these papers the authors 
fit their data to a universal scaling function which depends explicitly on li [18]. They argue that as the n alkanc 
increases from n = 11 to n — 14 the preferentially adsorbed component at the non-critical gas/liquid mixture surface 
changes from the n alkane to methyl formate and that this implies an effective surface field hi which is small and which 
changes sign as n is increased. More specifically they associated hi with the difference between the gas-liquid surface 
energies of the n alkane and methyl formate in the mixture. In practice this was taken to be the difference between 
gas-liquid surface tensions of pure n alkanc and pure methyl formate [18]. If hi determines whether a particular 
component is preferentially adsorbed at the interface it is certainly feasible that one can have hi = 0, at least for a 
certain range of temperatures. However, it is still not clear that a given n alkane-methyl formate mixture corresponds 
to a system undergoing an ordinary transition. An earlier reanalysis of surface tension data for isobutyric-acid-water 
mixtures suggests that the scaling properties are those of a system with a weak surface field, interpreted as the field 
that determines the preferential adsorption of one component at the non-critical interface [19]. The author also makes 
the important observation that similar weak-field characteristics should be found at the non-critical interfaces of any 
binary liquid mixture undergoing a wetting transition close to the critical endpoint [19,20]. For a one-component fluid 
adsorbed at a solid substrate this situation corresponds to the wetting transition being close to the gas-liquid critical 
point, a scenario which is expected if the wall-fluid attraction is weak. Crossover between strong and weak-field critical 
adsorption has been observed in very recent X ray diffraction studies of sublattice order in epitaxial FeCo films grown 
on an MgO substrate [21]; FeCo undergoes a continuous order-disorder transition from the B2 phase (CsCl structure) 
to the A2 phase (bec structure) and the substrate favors (weak) critical adsorption of the B2 phase for T >T C . 

Our paper is arranged as follows: in Sec. II we summarize results for the weak surface field behaviour of Ising 
systems in the semi-infinite lattice and in film geometry The magnetization profiles obtained for the latter are 
directly relevant to our study of a confined fluid. Sec. Ill reports the results of Monte Carlo simulations of the density 
profiles and adsorption of the truncated Lennard Jones fluid confined between two planar walls that exert a standard 
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(10-4) wall-fluid potential. Results are presented at the bulk critical point, and for some states on the critical isochorc, 
for various values of a parameter / that multiplies the (10-4) potential. By varying / we tune the strength of wall-fluid 
interactions relative to fluid-fluid interactions. / plays a role similar to that of hi for the Ising system and for certain 
values of this parameter we find the system is 'neutral'. This is confirmed by studies of capillary condensation below 
T c . In Sec. IV we describe the results of density functional calculations of the density profiles for the same model 
system as studied by simulation. By comparing the two sets of results we can assess the role of true surface criticality 
since the density functional approach is a mean-field one. We conclude in Sec. V with a discussion and summary. 

II. SUMMARY OF RESULTS FOR ISING SYSTEMS 
A. Semi-infinite lattice 

We consider an Ising magnet on a semi-infinite lattice with nearest-neighbour pairs of spins interacting (in bulk) 
with the exchange coupling constant J. The influence of the surface on the system is usually taken into account 
by introducing a modified coupling between spins in the surface layer, J\, and an (external) surface magnetic field 
Hi = J hi imposed on these boundary spins [1,2]. For vanishing bulk magnetic field, H = Jh = 0, the Hamiltonian 
of the model reads 

H = -J ^ a i a i ~ J i X! UiU i ~ Hl ^ Gi ( 2 - 1 ) 

<ij>ev <ij>esv iesv 

where V denotes the volume of the lattice and 5V the surface layer. In the absence of the surface field the leading 
critical behavior of various surface quantities depends on the ratio J/Ji and is described by three different surface 
universality classes [1,2]. At the bulk critical point, H = 0, T = T c , and at the special value Ji c ~ 1.50J the tendency 
to order near the surface is unchanged with respect to the bulk. This situation corresponds to the surface universality 
class of the "special transition". The condition Ji < J ic represents the universality class of the "ordinary transition" 
where the surface has a reduced tendency to order but becomes passively ordered at the bulk phase transition. For 
J\ > Jic spontaneous symmetry breaking also occurs above the critical point and on reducing T the bulk orders 
in the presence of an already ordered surface. This situation corresponds to the "extraordinary transition". In the 
renormalization-group analysis Ji c is an unstable fixed point, whereas for any starting value Ji < Ji c the surface 
coupling is driven to the stable fixed point Ji = by successive renormalization-group transformations. On the other 
hand, for Ji > Ji c the coupling is driven to the stable fixed point Ji = oo. In the framework of continuum field 
theory, appropriate to the universality class of the Ising model, the standard </> 4 Ginzburg-Landau model is augmented 
to include the surface by adding the surface contribution: 

HM) = \c<t? -hi4>, (2.2) 

where the parameter c can be related to the surface enhancement of the spin-spin coupling constant in the lattice 
models [2]. The leading critical behavior of a semi- infinite system with an ordinary or an extraordinary surface 
transition is described by stable renormalization-group fixed-point values c = +oo and c = — oo, respectively. The 
fixed point for the special transition is at c = c*. So far we have not included the effect of the surface field hi on 
the critical system. The presence of the surface field explicitly breaks the order-parameter (OP) symmetry at the 
surface inducing ordering in the surface layer for T > T c . It has been argued [3,4] that the extraordinary transition, 
i.e. hi — 0, Ji > Jic (c < c*), is equivalent to the case with hi ^ but with de-enhanced surface couplings: Ji < Ji c 
(c > c*). The latter case corresponds to the so called "normal" transition described by the fixed point hi = oo [4]. 
The term 'normal' is used because the situation with hi ^ and Ji < Ji c (c > c*) should be generic for simple 
fluids or fluid binary liquid mixtures bounded by a substrate or wall. Most substrates will attract molecules in a one 
component fluid giving rise to positive adsorption. Similarly for a binary mixture the substrate will favor one species 
of the fluid giving rise to preferential adsorption of that species. Both situations correspond to hi ^ 0. Moreover 
for fluids near walls the number of nearest-neighbor fluid 'bonds' is reduced with respect to the bulk which implies 
Ji < Jic or c > c* . 

A topic of considerable recent interest concerns the crossover region between the limiting values hi — oo and hi = 0. 
At the normal transition the OP (magnetization) profile m(z) decays monotonically from mi = 1 (all spins aligned) 
at the surface z — to the bulk equilibrium value which is zero for T > T c , the bulk critical temperature and is 
nonzero for T < T c . At the critical point the decay is described by a universal power law m ~ z~^l h ', where (3 and v 
are the usual bulk exponents, for distances large compared to all length scales. In d = 3 the Ising model exponents 
give f}/v = 0.518(7) [22]. By contrast, for hi = and T > T c , at the ordinary transition, the magnetization profile 
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is zero at any distance from the wall since the OP symmetry is not broken through the whole system [23]. It is of 
considerable interest to enquire about the profile for small but non-zero hi, i.e. in the crossover regime between the 
ordinary and normal transitions. Mean-field treatments do not predict any interesting behaviour of the magnetization 
profile; m(z) starts from a non-zero value mi and then decays monotonically to zero [1,3,24,25]. However, as has been 
shown [5] recently, for small hi close to the ordinary transition the critical fluctuations may lead to non-monotonic 
profiles. Near the surface m(z) may increase steeply with z to values much larger than mi. Near bulk criticality the 
range over which this rise occurs depends on hi through the length scale induced by this field. In an Ising system 
near the ordinary transition, i.e. for J <C Ji c , there is a single relevant surface scaling field hi [2] and the length scale 

induced by this field is l\ ~ h\ , where Aff d is the surface gap exponent [26]. At the bulk critical point the 

magnetization is given by 

m(z) ~ h lZ K for z < h (2.3) 

with k — (A° rd — (i)jv [5,7]. In the d = 3 Ising model A^ rd > j3 and n ~ 0.21 so m(z) increases with z near the 
surface. For z ~ l\ the profile has a maximum, and only for z > l\ does crossover to the 'normal' decay ~ z~^l v 
take place. For weak surface fields l\ can be large ( mesoscopic ) and thus the increase of the magnetization profile 
can occur over a mesoscopic distance. In turn, this may influence the behavior of the critical adsorption slightly 
above T c [6,17,18]. The short-distance behavior of the magnetization profile, Eq. (2.3), follows from the critical-point 
scaling of m(z) and the (continuity) assumption that m(z —> 0) ~ mi. Since at the ordinary transition the surface 
layer remains paramagnetic at T = T c it should respond linearly to a weak surface field hi [3]. The same should hold 
in the immediate neighbourhood of a surface, thus m(z) ~ mi ~ hi for z — > 0, hi — > 0. Assuming linear response 
for the scaling function leads to (2.3) [5]. A further heuristic explanation of the short-distance growth of the OP 
at the ordinary transition was given in Ref. [5,7] on the basis of the behaviour of correlations in the near-surface 
region. Using scaling analysis and arguments taken from short-distance expansions it was argued that there exists an 
effective parallel correlation length £,\\{z) which marks the crossover of the parallel correlation function from "bulk" 
to "surface" behavior. At T = T c the correlations decay algebraically. The decay in planes parallel to the surface is 
described by r ~ d+2 ~ Vs , where r is the distance within the plane parallel to the surface and the exponent r) s depends 
on the distance z from the surface. For z < l\, in the near-surface region, r) s assumes its bulk value ij (slow decay) 
if r <C z whereas for r ;§> z, r] s — rj^~ d > r\ (fast decay), so we can infer the effective range of correlations £|| ~ z. 
Now assume a weak surface field is imposed at bulk criticality on the surface with J\ <C J\ c - The magnetization 
induced in response at a distance z from the surface should be proportional to hi, to the correlation function in the 
perpendicular direction describing the decay of the order, < m(0)m(z) >^ z -<i+2--q ± ^ an( j ^ Q ^ ne corre i a t ec { area in 
the plane parallel to the surface, which can be influenced by a single surface spin. Together these terms give 

m(z) ~ /ii < m(0)m(z) > - hiz 1 -^ for z < h (2.4) 

Scaling relations [7] give in turn 1 — ■q°^ d = k and we recover (2.3). Many of the above predictions were confirmed by 
Monte-Carlo simulations for the d — 3 Ising model [8] and by detailed renormalization group studies of a cf> 4 model [7]. 
Exact results [10-12] and Monte Carlo simulations [9] for the d — 2 Ising model yield magnetization profiles that also 
increase near the surface for the case of a weak surface field. However, in d = 2 the simple power law (2.3) is modified 
by a logarithm, i.e. m(z) <~ h\z K \og(h\z), where k = 3/8. Note that in mean- field theory Af d = (3 = 1/2 so that 
k = and the magnetisation profile decays monotonically as stated above. 

B. Ising films 

In this subsection we consider film geometry, i.e. a lattice of finite width D in the z direction but infinite in the 
other directions. The critical properties of such a confined system should be particularly sensitive to the value of the 
(identical) surface fields when the length scale l\ becomes comparable to or larger than the width D. Recent studies 
of critical d = 2 Ising films of finite width D subject to identical surface fields h\ = ho using exact transfer-matrix 
diagonalization and density-matrix renormalization-group (DMRG) methods fully confirmed this expectation [13]. 
For weak surface fields such that l\ < D, magnetization profiles at T — T c , h — have two symmetric maxima shifted 
away from the walls and lying at z ~ l\ and z ~ D — l\ (see Fig.l). In this regime of hi the surface magnetization 
mi = mo decreases fairly rapidly as the value of hi decreases, whereas the value of the magnetization at the midpoint 
of the film, m(D/2), changes only very slightly. The two separate maxima merge into one, broad maximum located 
at the middle of the film for smaller hi where li ~ D. For even smaller hi, li > D and the profiles are nearly flat with 
the magnetization much lower near the walls than in the central part of the system. In this regime the magnetization 
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both at the surface and near the midpoint, in fact through the whole profile, decreases as ~ h\ h\h\ (the logarithmic 
factor is specific to the d = 2 Ising model). This result agrees with the physical interpretation of l\ as a length scale 
up to which the OP profile responds linearly (up to the logarithmic factor) to the surface field. Note that the change 
in the shape of the magnetization profile from one with two maxima to one with a single maximum at the center of a 
film can be achieved either by decreasing h\ at fixed D or by decreasing D at fixed h\. This change is accompanied 
by a rich variation of the solvation force [13]. The latter is an important thermodynamic quantity characterizing a 
confined system [27]. For a fluid the solvation force is the excess pressure over the bulk arising from confinement by 
walls; it can be measured by the surface force apparatus. For the Ising system the solvation force = —{df*/dD)T,hi,h, 
where f *(D) is the finite-size contribution to the free energy The critical-point scaling function of the solvation force, 
i.e. at T = T c and h = 0, has a pronounced maximum near l\ — D (see Fig. 6 of Ref. [13]). 



III. SIMULATION STUDIES OF THE LENNARD-JONES FLUID IN A SLIT PORE 

In this section we describe the results of a Monte Carlo study of a three dimensional system, namely the truncated 
Lennard-Jones fluid, confined between two planar walls separated by a distance D. By decreasing the strength / of 
the wall-fluid interactions while keeping the fluid intcrparticlc potential constant we examine the effects of decreased 
wall-fluid attraction on the shapes of the density profiles and on the form of the adsorption in the critical fluid. 
Decreasing / should, in some sense, correspond to reducing h\ in the Ising system discussed in the previous section. 



A. Computational details 



We have performed Monte Carlo simulations of the Lcnnard-Jones (LJ) fluid having interparticle interactions of 
the form: 



U L .,(r) = 4e 



(3.1) 



Here e measures the well depth of the potential, while a sets the length scale. The potential was truncated at 
a radius r c = 2.5a and left unshiftcd. No corrections were applied to account for effects of the truncation. The 
simulations were performed within the grand canonical (constant-^ FT) ensemble [28] , permitting fluctuations in the 
total particle number N. In order to mimic a slit pore geometry, a cuboidal simulation cell of dimensions L x x L y x D 
(with L x = L y = L) was employed. Structureless planar walls were imposed in the planes z = and z = D, while 
periodic boundary conditions were applied to the cell boundaries in the x and y directions parallel to the walls. Fluid 
particles were assumed to interact with a single planar wall via a long range potential having the form: 



U w (z) = 4ef 



icir-cf)' 



(3.2) 



where / is a parameter that tunes the strength of the fluid- wall interactions relative to those of the fluid interparticle 
interactions. Decreasing / simply reduces the depth of the minimum. The total wall-fluid potential is U w {z) + U W (D — 
z). We note that U w (z) models a wall comprising a single plane of LJ particles [29]. No potential truncation was 
applied to U w (z) since the wall potentials decay considerably less rapidly with increasing separation than the LJ 
interparticle potential of eq (3.1). Further details of the simulation procedure are given in an earlier study [30]. In 
the course of the simulations the one-dimensional density profile was accumulated: 



P (Z) EE f 
JO 



p(r)dxdy 



(3.3) 



representing the configurationally averaged local number density at a given distance z. Measured forms of this profile 
for a range of wall strengths / at the liquid-gas critical point and along the critical isochore are presented in the 
following subsection. 



B. Density profiles at the critical point 

The bulk critical point parameters of the truncated LJ fluid have previously been determined from an accurate 
finite-size scaling study [31]. They arc (in standard reduced LJ units [28]), T c = 1.1876(3), /J, c /k B T c = -2.778(2). 
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Using these parameter values we have obtained the density profiles p(z) for a range of wall strengths /. The results 
for a system of size L — l5a,D — 40cr are shown in Fig. 2(a). A magnified plot of the region around the critical 
density p c — 0.3197(4) (horizontal line) is shown in Fig. 2(b). Examination of the profiles of Fig. 2 reveals three 
qualitatively distinct regimes of behaviour as the wall strength / is varied. We discuss them in turn. For the smallest 
values of / studied, the density near the wall lies well below its critical value and the profile exhibits little structure. 
Moving away from the wall, however, the density increases, reaching a maximum in the slit middle, although it at no 
stage attains the bulk critical value. The overall shape of the profile is convex upwards with respect to the critical 
density. For the largest values of / studied, the density near the wall is much higher than the critical density and 
large packing effects are evident in the profile. As the distance from the wall increases, the density falls progressively, 
although it does not reach the bulk value. The profile is convex downwards with respect to the critical density. For 
intermediate values of /, the packing-induced density oscillations near the wall span the critical density, while further 
away from the wall the magnitude of the profile curvature is generally less than in the large or small / limits. A 
magnified view of the region around the critical density is given in Fig. 2(b). Strikingly there exist some / values in 
this regime for which p(z) exceeds p c , but the profile is concave upwards with respect to the critical density. Thus 
p{z) exceeds p c and increases with increasing distance from the wall. As / increases, however, one sees (Fig. 2(c)) 
the formation of a structure with two smooth maxima in the profile, the maxima being some distance from the walls. 
On further increasing /, the maxima become more pronounced and their positions move closer to the walls. Clearly 
these results show very similar features to the magnetisation profiles shown in Fig.l. ft is also instructive to observe 
the effect of changing the slit width on the double maxima profile of Fig. 2(c). The comparison with a profile for a 
slit of width D = 20a is made in Fig. 3, from which one sees that the double maxima profile is replaced by one with 
a single maximum in the narrower slit. Once again this is similar to what is found for the Ising system (Sec. 11). 

C. Properties of the Neutral Wall System 

The results of the last subsection suggest that for certain choices of the strength parameter / the confined Lennard- 
Jones fluid is behaving similarly to an Ising film with hi « 0. Density profiles for / < 0.653 are of particular interest 
since it is in this regime of / that there is a crossover from a maximum to a minimum in the profile at the midpoint. 

One of the signatures of the ordinary transition (at a single surface) is that the magnetisation profile remains 
constant, m; = 0, 1 = 1,2, ■•■ for all T > T c and h = 0. The analogous situation for a confined fluid should be 
p(z) = p c , the critical density, except for distances close to the wall where packing constraints induce oscillations 
in the profile, for states on the critical isochore: pbuik(n,T) = p c ,T > T c . This implies, in turn, the adsorption 

r = J Q (p(z) — p c )dz should be constant (and close to zero) along the critical isotherm. 

We varied / in an attempt to find a 'neutral' wall where such properties were achieved. Figs 4(a) and 4(b) display 
the density profiles on the critical isochore for / = 0.644, D = 40c and 20ct, respectively. In each case the density 
profile is almost flat throughout the slit except for near the walls. Fig. 5 shows a plot of the adsorption To 2 versus 
t = (T — T c )/T c for D = 20a and a selection of values of /. We observe that T is temperature independent for the 
'neutral' wall with / = 0.644. Increasing / slightly leads to positive adsorption near T c whereas decreasing / below 
0.644 leads to more pronounced negative adsorption which is T dependent. 

Note that the simulations were performed at chemical potentials p measured along the critical isochore of the bulk 
fluid. These were determined, using periodic boundary conditions, in Ref. [30]. 

Further evidence that / = 0.644 corresponds to a 'neutral' wall is furnished by considering the confined fluid below 
T c . For strongly attractive walls, favouring the liquid phase, we expect liquid-gas coexistence to occur at a chemical 
potential p which lies below the value p cx for bulk coexistence at the same temperature [27] whereas for purely 
repulsive walls, favouring gas, coexistence should occur at a higher value than in bulk. We determined the line of 
capillary coexistence , where liquid and gas coexist, in the slit with D — 20a. For —0.1 < r < and / = 0.644 we 
found this line was hardly shifted from the bulk coexistence line - see Fig. 6. The greatest deviation, at r = —0.1, was 
about 0.1%, with condensation in the slit ocurring for a (slightly) higher value of p, implying the walls with / = 0.644 
favor (slightly) the gas phase, i.e. the contact angle is (slightly) larger than n/2 at this temperature [27]. The inset to 
Fig. 6 displays the difference between the chemical potential p, for capillary coexistence, and the corresponding value 
p cx , for bulk coexistence, for three values of /. One observes that for / = 0.71, p — p cx remains substantially more 
negative than the result for / = 0.644 whereas for / = 0.59, corresponding to a weaker wall attraction, p — p cx is 
positive at all temperatures in the range considered. These results suggest that a value of / slightly lower than 0.644 
would yield capillary coexistence very close to that in bulk, i.e. would correspond in magnetic language to a surface 
field hi = as regards the phase behaviour below T c . 
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IV. DENSITY-FUNCTIONAL THEORY STUDIES OF THE LENNARD-JONES FLUID IN A SLIT PORE 



Having established that computer simulations of the critical Lennard-Jones fluid do find features in the density 
profile that mimic those found for the magnetisation profile in systems with weak surface fields hi, we now enquire 
whether any of these features are also found in a mean-field treatment of the same confined fluid. The oscillations 
in the profiles shown in Figs 2-4, arise from packing effects near the walls and these should certainly be captured by 
a suitable density functional theory (DFT) approach. Recall, however, the double maxima, separated by a shallow 
minimum, is a feature which reflects true surface criticality in the case of a magnetic system (Sec. II). Does the 
presence of a similar feature in the simulation results for fluids also signal some manifestation of surface criticality 
or is the feature somehow reflecting the particular form of the wall-fluid potentials? By carrying out mean-field 
calculations we attempt to address these questions. 

We consider the same system of particles as in Sec. Ill, i.e. the Lennard-Jones fluid with the pair potential given 
by (3.1), but not truncated, and apply the DFT [32] to determine the density profile at the bulk critical point. All 
the functionals implemented here are of mean- field type; bulk critical exponents would take their mean-field values. 
We assume that the fluid is confined between two identical parallel walls located at z = and z = D, and infinite in 
the x and y directions. As previously, the system is in contact with a reservoir of fluid at temperature T c and critical 
chemical potential p c . The grand potential of the inhomogeneous system is a functional of the one-body density 
distribution p(r) [32]: 

Q\p] = T[p] - J dr(p - V(r))p(r), (4.1) 

where V(r) = V(z) is the total wall-fluid external potential and the equilibrium density profile p(r) = p(z) corresponds 
to the minimum of The simplest form for 0[p] which should be appropriate near criticality is based on the square 
gradient approximation to the intrinsic Helmholtz free energy functional T[p] with the wall-fluid contribution modeled 
by a term <J> S which depends only on the fluid density at contact i.e. $ s = |(/9 2 (0) + p 2 {D)) — e w (p(0) + p{D)), where 
c and e w are constants characterizing the effects of the wall. Although functionals of this type cannot incorporate 
short-ranged correlations arising from packing and hence cannot acount for the oscillatory behavior of the density 
profile near the walls [27], they were succesfully employed in studies of finite size effects on the critical adsorption in 
large pores, as this phenomenon is dominated by the slow decay of the profile far from the walls [33,30]. However, 
the density profiles obtained from this approximation are always monotonic; depending on the value of c and e w , p{z) 
is either monotonically increasing or monotonically decreasing. This observation follows straightforwardly from the 
graphical construction for the solution of the Euler-Lagrange (differential) equation [33,34] and remains valid for the 
local functional approach employed by Borjan and Upton [35] in their studies of a critical confined system. Thus it 
is necessary to consider other approximations for the intrinsic free energy functional [32,36], in which T[p] is divided 
into contributions to the free energy functional arising from repulsive forces, T r \p\, and from the attractive part of 
the pairwise potential; the latter is treated in mean-field or random-phase approximation. Using hard-spheres (HS) 
as the reference fluid one constructs a Van der Waals approximation: 

T\p\ = T HS [p\ + \ J dri J dr 2 p(ri)p{Y 2 )<p att (n 2 ), (4.2) 

where ri 2 = |ri — r 2 |. The attractive part of the pairwise potential is not specified uniquely. Here we employ the 
prescription: 

4>att = U LJ (r ) r < r 

= U L .,(r) r > r 0) (4.3) 

where r = 2 1 ^ 6 a is the minimum of the LJ potential and Ulj{to) = — e. In this approximation the Helmholtz free 
energy density of a uniform fluid has the form: 

M = f HS (p) - Y a > ( 4 - 4 ) 

with a = — J dr<j> at t(r) > 0. For simplicity we set the hard sphere diameter d equal to a. We describe results of 
calculations using both the weighted density approximation (WDA) and the local density approximation (LDA) for 

T HS [p\m. 
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A. Results from WDA 



In order to incorporate correctly the main features of the short-range structure of a fluid near the walls, we employ 
a simple non-local approximation for Ths [p] based on the weighted-density approximation [37-39] : 

Fhs[p\ = ?id[p] + j drp(r)Vex(p(r)), (4.5) 

where the ideal gas contribution is 

T ld [p\ = J drp(r)[ln(A 3 p(r)) - 1]. (4.6) 

A is the thermal de Broglie wavelength. vp ex is the excess, over ideal, of the free energy per particle of the uniform 
fluid. In (4.5) this quantity is calculated at some weighted density p, where p is a non-local functional of the true 
density p obtained by taking the weighted average of p(r) over a local volume reflecting the range of interatomic 
forces: 

p(r) = J dv'p{v')w{\v - r'|). (4.7) 

The precise definition of w depends on the particular version of WDA; the most sophisticated versions introduce 
weight (s) so that the second functional derivative of Ths[p] reproduces the two-particle direct correlation function 
c( 2 )(r) of a uniform hard sphere fluid [36]. In this paper we assume the simplest, density independent step function 
form for w(r), i.e. w(r) — 3/(47rd 3 )0(<i — r), where Q(x) denotes the Heaviside step function. This choice yields 
reasonable results for the density profiles near a hard wall [38]. The excess free energy per atom is given by the 
Carnahan-Starling expression: 

IMP) = k B T V ^~^ 2 \ (4.8) 

where r\ = \itpd? is the packing fraction. Minimizing Q[p] in (4.1) with the external potential V(z) yields an integral 
equation for the density profile of the fluid: 



p{z) = po exp{-(l/fc s T)hMp(*)) + J dv' p{z')^ ex {p(z'))w{\v - r'|) 

+ J dr'p(z')ci> a u(\r-r'\) + V(z)}}, (4.9) 

where /?o = A 3 exp(p/ksT). 

We solve this equation numerically for two choices of the single wall potential U w {z): (a) a long range potential 
with the same form used in simulations (3.2) and (b) a hard wall potential with an attractive short-ranged tail: 

U ws (z) = 00 z < 

= -efexp(-zfa) z>0. (4.10) 

Once again the strength parameter / determines the depth of the potential well. The integral equation (4.9) was 
solved using a minimization procedure [40] which, in principle, is always convergent; the commonly used iteration is 
much less efficient near bulk criticality where profiles decay very slowly. We minimimize the function 

N 

u(pi, ■ ■ ■ , Pn) = ^2(Pi - Pi? ', (4-11) 

i=i 

where pi — p(i) are densities calculated at N points between z = 1 and z = D and are treated as independent 
variables, pi denotes the right-hand side of the integral equation (4.9) calculated at z = Zf, pi are functions of pi, i.e. 
Pi = PiiPii ■ ■ ■ ! Pn)- The minimum u « corresponds to the solution of (4.9). This minimum is determined using the 
conjugate gradient method. Calculations were performed for the slit pore of the same width as in the simulations, 
D = 40(7, at the critical point of the uniform fluid. The mean field free energy density (4.4) gives, for our model, 
ksTc/e = 1.00617(2) and p c a 3 — 0.24912(4). To allow for the expected oscillatory behavior of the profile near the 
walls we used a very fine mesh, N — 400. As a starting density profile we usually took p(z) = p c . 
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In Fig. 7(a) we show density profiles obtained for the long-ranged wall-fluid potential (3.2) for various strengths / 
between and 1. The sequence of profiles is very similar to that obtained in simulations (see Fig. 2(a)), i.e. there 
are three different types of profile shape depending on the value of /. For the smallest values of / studied ( between 
and about 0.55), the profiles are convex upwards with respect to p c and exhibit little structure, whereas for the 
largest values of / studied (between about 0.7 and 1) the profiles are convex downwards with p(z) > p c and show very 
pronounced oscillations near the walls. At intermediate values, i.e. as / increases above w 0.55 the profiles become 
flatter in the central part of the slit. They remain concave upwards with respect to the critical density but move as 
a whole towards greater values so that except near the walls, p(z) > p c - see Fig. 7(b). A difference between results 
from the WDA and those from simulations occurs when the two smooth maxima, separated by a shallow minimum, 
structure forms in the profile. In the simulations these new maxima are more pronounced than in WDA and are 
located somewhat closer to the walls. It is instructive to compare the height and position of the new maxima with 
those of the third peak in the density profile, arising from packing effects (layering) near the walls. In the WDA this 
third peak is located at z w 3<r. For / < 0.55 the peak is not pronounced indicating weak packing effects in the 
weakly attractive wall potential. As / is increased the third peak becomes more pronounced and the value of the 
density at this peak increases rapidly. For / < 0.63 this density lies below the density at the midpoint, p m id, and 
below the density at the new maxima (where these have developed). However for / in the range 0.63 < / < 0.64 the 
density in the region of the third peak grows more rapidly than in the central portion of the density profile and for 
/ > 0.64 the density at the third peak is higher than at the new maxima and p m id- The positions of the new maxima 
move closer to the walls as / increases and for / > 0.68 these maxima cannot be discerned. Note that the range of 
values of / corresponding to the three qualitatively different regimes of behavior of the profile are very close to those 
found in simulations. 

In Fig. 8 we plot the adsorption Fa 2 versus t, for D = 40cr, calculated along the critical isochore p = p c , for several 
values of /.These results should be compared with those from simulation in Fig. 5 noting that the latter correspond 
to a much wider range of the strength parameter /. We should also note the difference in vertical scale between these 
two figures; the adsorption is very weak within the WDA for the narrow range of / displayed in Fig. 8. On the basis 
of this figure we might conclude that a value of / in the range 0.561 < / < 0.564 corresponds to a 'neutral' wall. In 
this range T changes from monotonically increasing with r to exhibiting a shallow minimum. There is some value for 
which T(r) is flat and negative, c.f. Fig. 5. 

Figure 9(a) presents a selection of the profiles calculated for the short-ranged wall-fluid potential (4.10) with various 
strengths / between 1 and 0. There are generally fewer oscillations near the walls and the contact values p(0) = p{D) 
are large for large values of / but away from the walls the shape of the profile changes in a similar way as for the 
long-ranged U w (z) as / is increased. The two smooth maxima structure forms at much larger values of / however, 
i.e., at about 0.8 and then shifts rapidly towards the walls, merging with oscillations - see Fig. 9(b). Such behaviour 
might reflect the narrower attractive well in the wall-fluid potential. One requires a larger / to acquire the same 
strength of attraction. In the range of / we studied there are no profiles characteristic of strongly adsorbing walls 
with extreme packing effects. 



where fnsip) is the Hclmholtz free energy density of a uniform hard-sphere fluid with density p [32,36]. The LDA 
does not incorporate the short-ranged correlations that characterize the structure of dense liquids and hence it cannot 
describe oscillatory profiles arising from packing effects. Attractive interparticle forces are still treated by the mean- 
field approximation given in (4.2) and the free energy density of the uniform fluid is still given by (4.4) so that the 
critical point remains the same as that given in the previous subsection. 

Critical point profiles for the slit of width D = 40a obtained in the LDA for the long-ranged wall-fluid potential 
(3.2) are shown in Fig. 10(a). The wall strengths / are the same as in Figs. 9. Despite the very crude approximation, 
some weak structure remains in the profiles. The sequence of profiles is reminiscent in shape to those from Fig. 9, 
i.e. for WDA with the short-ranged wall-fluid potential. For / < 0.8 the shapes are convex upwards with respect to 
the critical density and, except for the first peak near the walls, p(z) lies below p c . Profiles for / > 0.8 are similar in 
shape to these obtained in simulation and in WDA for intermediate values of / - see Fig. 10(b). As was the case for 
WDA with the short-ranged wall-fluid potential, the characteristic two smooth maxima structure is formed for larger 
values of / than in simulations. 



B. Results from LDA 



In the local density approximation (LDA) the free energy functional of the hard-spheres has a form: 




(4.12) 
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V. DISCUSSION 



In this paper we have performed extensive Monte Carlo simulations and density functional calculations of the density 
profile and Gibbs adsorption T for a Lennard- Jones fluid confined between two planar walls under conditions where the 
reservoir is at, or close to, bulk criticality. We observe that increasing the parameter / that multiplies the (10-4) wall- 
fluid potential alters the profile dramatically, from one typical for a purely repulsive wall (T is strongly negative) to 
one typical of a strongly attractive wall (T is strongly positive), for the critical fluid at fixed wall separation D. In the 
simulations we find that for / « 0.644 the confined LJ fluid behaves in a similar fashion to an Ising film with vanishing 
surface field hi = 0, i.e. to what is observed at the ordinary transition. This 'neutral' wall scenario is characterized 
by a fluid density profile which, away from the walls where oscillations arise, is almost constant throughout the slit 
and which gives rise to a T that is constant along the critical isochore (see Figs 4 and 5). The same value / « 0.644 
gives rise to a line of capillary coexistence for T < T c that differs very little from the bulk coexistence line (see Fig. 6), 
supporting the contention that this value corresponds to a 'neutral' wall. Moreover, we observe that as / decreases 
from 1 the density profile of the critical LJ fluid, away from the walls, exhibits features reminiscent of those found for 
the magnetisation profile of the critical Ising film in the crossover regime between the 'normal' and ordinary surface 
universality classes which is achieved by reducing the surface field hi towards zero. In particular for 0.7 > / > 0.64 
and D = 40a we find profiles with two smooth maxima, located away from the walls; compare Figs 2b and 1. As / 
is reduced further, or as the slit width D is reduced at fixed /, these maxima merge into a single maximum located 
at the midpoint D/2. In the Ising film a magnetisation profile with two smooth maxima is characteristic of a system 

subject to weak surface fields hi = ho] the profile increases for distances z < h ~ h x 1 < D/2, before decaying 
to the value at the midpoint (see Sec. II). For weaker hi (or smaller D) a single maximum forms at the midpoint 
when li ~ D/2. 

In the light of the simulation results it is tempting to argue that the near critical LJ fluid does exhibit several 
features in common with the Ising system. Consideration of the adsorption and of capillary condensation would 
suggest identifying hi = f — f neutral, with f neutral ~ 0.644. However, the fluid density profiles at the critical point 
clearly do not exhibit the symmetry m(z; —hi) = —m(z;hi) of the magnetisation profile of the Ising system - even 
when some attempt is made at coarse graining the oscillations near the walls and subtracting a constant density. Thus 
it is not obvious what is the best choice of order parameter profile for the fluid. The form of the confining potential 
obviously plays an important role and in order to investigate this further we performed a series of DFT calculations 
using a very simple WDA to incorporate short-ranged correlations. 

Our DFT results also point to the existence of a 'neutral' wall; the adsorption T(t) calculated along the critical 
isochore is almost fiat for / lying between 0.561 and 0.564 pertaining to the same (10-4) wall potential and wall 
separation D = 40cr as the simulations - see Fig. 8. Moreover, the DFT density profiles exhibit the characteristic 
two smooth maxima for an intermediate range of / observed in the simulations - see Fig. 7. What does this lead us 
to conclude about the role of critical fluctuations in determining the shape of the density profiles for 'weak' surface 
fields? Recall that unlike simulations DFT is a mean-field treatment which does not incorporate effects of fluctuations. 
The critical scaling arguments of Sec. II state categorically that the mean field order parameter profile for the Ising 
system subject to a surface field hi acting on the surface layer only should vary monotonically away from the wall, 
i.e. there should be no maxima. The mean-field value of the exponent that describes the short-distance behaviour of 
the profile, Eq. (2.3), is k = 0. That we find smooth maxima in the density profiles within our DFT treatment would 
seem to reflect the shape of the confining (wall-fluid) potential. Note that the precise form of wall-fluid potential 
is not crucial as we also find the maxima, albeit for larger values of /, for a short-ranged, exponentially decaying 
attractive potential - see Fig. 9. There are quantitative differences between DFT and simulation. The maxima are 
more pronounced, are located closer to the walls and appear to shift faster towards the midpoint as / is decreased in 
simulation and this might reflect the presence of fluctuation effects. The latter give rise to a slightly higher value of 

the exponent determining the length scale l\ <~ h 1 v ^ Al than in mean field, i.e. the d = 3 Ising value is 1.3 rather 
than 1. These observations suggest that one should re-analyse the arguments which lead to Eq. (2.3). Both arguments 
given in Sec. IIA are based on the assumption that the surface magnetisation mi and m(z ~ 0) respond linearly to 
a weak hi for T = T c . Whilst this is eminently plausible for the Ising system it is not clear how to translate to the 
case of a fluid near a wall. It is possible that the shape of the external potential has non-trivial repercussions for 
the short-distance expansion of the order parameter profile. However, one should also bear in mind that there is a 
further complication in attempting to relate crossover behaviour in fluids to that in systems with Ising symmetry 
The reduced order parameter symmetry of bulk fluids lead to scaling field mixing [41,42]. However, to our knowledge, 
there have to date been no considerations of the consequences of this reduced symmetry in the context of surface 
critical behaviour of fluids. Any scaling theory that seeks to encompass fluid and magnet surface critical behaviour 
within a unified description must address this issue. 

One way of doing so may be to generalize the role of the surface scaling fields c and hi (defined for Ising-like 
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systems) such that regardless of the model system, they correspond to equivalent paths in the surface phase diagram. 
Consider the case of an Ising model in zero bulk field at its bulk critical temperature. Within the surface phase 
diagram, the line hi = corresponds to a special symmetry line of the system because the magnitude of the surface 
order parameter is invariant with respect to a change of sign of hi. The scaling field c (which relates to missing 
neighbours and surface induced modifications of surface couplings), controls the temperature at which the surface 
orders relative to the bulk. The sign or magnitude of c does not determine which phase is preferred at the wall, thus 
the c field can be represented as a vector coincident with the symmetry line hi — (Fig. 11(a)). 

By contrast in a fluid, missing neighbours typically lead to a segregation of one phase to the wall, i.e. to breaking 
of the order parameter symmetry with respect to its bulk critical point value. In terms of the geometrical picture 
outlined above, the line of surface order parameter symmetry lies off-axis in the surface phase diagram and does not 
coincide with the c-field. In attempting to restore the symmetry, one might consider applying a non-zero external 
surface field hi that couples to the order parameter near the surface, and whose magnitude is chosen such as to return 
the value of the order parameter at the surface to its bulk critical value. Within this framework, the combined c and 
hi fields can be regarded as an effective field c directed along the symmetry line. Similarly we require an effective 
field hi (the scaling equivalent to the Ising model hi) which measures deviations from this line (Fig. 11(b)). Formally 
these ends can be attained by introducing 'mixed' surface fields having the form 

c = (c — c*) + aihi (5-1) 
hi = hi + a 2 (c - c*) , 

with ai and a 2 system specific constants (both zero in the Ising model context) whose values control the degree of 
field mixing. 

The starting point for exploring the consequences of field mixing for surface critical behaviour is an expression for 
the singular part of the surface free energy. In terms of the mixed fields this is 

F! ur « f 2 - a ff ur (ilf- A ,cf-^hif- Al ) , (5.2) 

where /j^ r denote scaling functions, while f and \X are (mixed) bulk scaling fields which in the Ising context reduce to 
the reduced temperature and reduced bulk field respectively [42]. The positive constants A and Ai are respectively 
the gap exponent and surface gap exponents, while <fi is an exponent controlling the crossover from the ordinary to 
the normal universality classes. 

We shall not attempt to explore the consequences of surface field mixing in any great depth here, although we do 
hope to report on the matter in future work. For the time being, we merely point out that one corollary of Eq.(5.2) 
is that-in a simple fluid-the scaling equivalent of the Ising model magnetisation profile m(z) is not the density 
profile p(z), but a profile formed by linearly combining the density and the energy density profiles. One should also 
consider the adsorption not on the critical isochore, but on an analytical extension of the coexistence diameter. These 
observations may prove important in exposing more fully the universality of surface critical phenomena for fluids and 
magnets. 

Finally, with regard to the applicability of the above considerations to realistic fluids, the question arises as to 
just what constitutes the surface field hi in a fluid. Clearly the wall-fluid potential in general couples to the order 
parameter in the vicinity of the surface and thus can be used to tune its value. However, as noted in Sec. I, a realistic 
potential such as the 10-4 potential used in this work will not allow one to restore the order parameter to its bulk 
critical value exactly at all distances from the wall. Short ranged packing effects will always be visible close to the 
wall. Nevertheless, in the context of the ordinary transition, the length scale on which packing effects occurs will 
be negligible compared to the length scale of the critical phenomena (set by h). Hence one might speculate that 
fluid-magnet universality should be evident provided one views the system on sufficiently large length scales, i.e. such 
that the value of a suitably coarse-grained order parameter profile matches its bulk value near the wall. 
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FIG. 1. Magnetization profiles mj of the d — 2 Ising film of width D = 200 lattice spacings for T = T c and zero bulk field and 
several values of the surface field h 1 =h D : from top to bottom profile: hi = 10, 0.5, 0.4, 0.3, 0.2, 0.14, 0.1, 0.07, 0.04, 0.01, 0.001. 
These results were obtained using exact transfer-matrix diagonalization (adapted from Ref. 13). 
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FIG. 2. (a) The measured critical point density profiles p(z) (in units of cr) for a LJ system of size L — 15a and width 
D = 40a. Wall strengths / from top to bottom: / =1.0, 0.9, 0.8, 0.7, 0.68, 0.67, 0.66, 0.653, 0.64, 0.63, 0.62, 0.5, 0.4, 0.3. The 
horizontal dashed line denotes the bulk critical density, p c = 0.3197(4) (b) Magnified version of (a) showing the region close to 
the critical density. Wall strengths / from top to bottom:/ = 0.68, 0.67, 0.66, 0.653, 0.64, 0.63, 0.62. (c) Magnified version of 
the profile for / = 0.653 showing the two smooth maxima with a shallow minimum in the centre of the slit. 




FIG. 3. Critical point density profile p(z) for / = 0.653 for slits of width D — 20a and D — 40a. In both cases L = 15a. 
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FIG. 4. Density profiles on the critical isochore for / = 0.644, the 'neutral' wall, (a) L — 15a and width D — 40a (b) 
L — 15a and width D = 20a. Note how insensitive the profiles are to the value of T* 
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FIG. 5. The adsorption T, in units of a 2 , calculated along the critical isochore as a function of reduced temperature 
r = (T — T c )/T c for various values of /. Results correspond to L = lba and width D = 20a. Note that for the 'neutral' 
wall / = 0.644, T(r) is constant. 
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FIG. 6. The chemical potential /i (in units of fcgT) versus reduced temperature r for capillary condensation in a system of 
size L = 15a and width D = 20<r for wall strength / = 0.644. The thick solid line denotes the bulk coexistence line ficx- The 
inset displays the difference A/i = fj, — ji cx between the coexistence lines for three values of /. For / = 0.59, /i — jj, cx is positive, 
i.e. there is capillary evaporation. 
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FIG. 7. Critical point density profiles p{z) (in units of u 3 ) for a slit of width D = 40ct and the long-ranged wall-fluid 
potential U w (z) (3.2) obtained using the WDA for several values of the wall strength /: (a) from top to bottom: / =1.0, 0.9, 
0.8, 0.7, 0.65, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1. (b) Magnified version for /= 0.66, 0.65, 0.64, 0.63, 0.62, 0.6, 0.55 and 0.5. Note 
that for 0.62 < / < 0.68 two smooth maxima, separated by a shallow minimum at the midpoint, arise at distances away from 
the walls, outside the region where the oscillations due to packing occur. The horizontal dashed line denotes the bulk critical 
density, p c = 0.24912(9). 
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FIG. 8. WDA results for the adsorption T, in units of a 2 , calculated for a slit of width D = 40a along the critical isochore 
as a function of reduced temperature t, for several values of /. The wall-fluid potential U w (z) (3.2) is long- ranged. Note that 
the adsorption is very small and only weakly dependent on temperature for the range of values of / shown here. 
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FIG. 9. Critical point density profiles p(z) (in units of <r 3 ) for a slit of width D = 40a and the short-ranged wall-fluid 
potential U wa (z) (4.10) obtained using the WDA for several values of the wall strength /: (a) from top to bottom: / =1.0, 
0.9, 0.85, 0.8, 0.7, 0.68, 0.66, 0.64, 0.62, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1. (b) Magnified version of (a) for the largest values of / : 1, 
0.9, 0.85, 0.8, 0.7, 0.68, 0.66. The horizontal dashed line denotes the bulk critical density, p c = 0.24912(9). 
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FIG. 10. Critical point density profiles p(z) (in units of a 3 ) for a slit of width D — 40a and the long-ranged wall-fluid 
potential U w {z) (3.2) obtained using the LDA for the same choice of wall strengths / as in Fig. 9(a). (b) Magnified version for 
/ = 1,0.9, 0.8. The horizontal dashed line denotes the bulk critical density, p c = 0.24912(9). 
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FIG. 11. The form of the surface scaling fields c and hi, for a system at its bulk critical point parameters, as described in 
the text, (a) A model of the Ising symmetry where there is no mixing, so c = c and hi = /ii(b) A fluid system. 
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